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PARAMETRIC STUDY OF A THERMOELECTROSTATIC GENERATOR 


FOR SPACE APPLICATIONS 
by Paul M. Margosian 
Lewis Research Center 


SUMMARY 

Examined in this report is the theoretical performance of a thermoelectro- 
static generator, which is a variable capacitance device capable of converting 
thermal energy directly to electrical energy. Two physically realizable ther- 
modynamic cycles were selected. Two dielectric materials, barium titanate and 
polyethylene terephthalate (Mylar), were selected for their potential ability 
to yield the best performance attainable with this device. Mathematical ex- 
pressions for internal energy and entropy in terms of electrical and thermo- 
dynamic parameters were derived. Mathematical approximations to the behavior 
of the dielectric constants of the materials were chosen. Equations describing 
the thermodynamic cycles were derived and solved numerically. The upper limits 
on performance were found to be the following: maximum thermal efficiency, 

1 to 4 percent; and maximum energy output, 0.5 to 1.5 joules per cubic centi- 
meter of dielectric per cycle. Realistic performance, dictated by material 
properties and limitations on the electric field that can be supported, is 
lower by at least a factor of five. 

For use in space an electrical generator should have minimum weight. In 
order to evaluate this parameter, the heat-transfer equations were written for 
a typical design, heated by the sun or by a nuclear reactor, and solved numer- 
ically. Combining the results of the thermodynamic and the heat-transfer cal- 
culations led to an estimated minimum specific weight of the order of 50 to 
150 pounds per kilowatt. 


INTRODUCTION 

Lightweight sources of electrical power are essential for the performance 
of almost any mission in space. Power requirements range from a few watts for 
communication with near-Earth satellites to several megawatts for a Mars mis- 
sion employing electric propulsion (ref. l) . 

An attractive method of producing electric power in space is the direct 
conversion of thermal energy, such as solar radiation or heat from a nuclear 
reactor, to electric energy. One such method is the thermoelectrostatic gener- 




Silicon 

controlled 

rectifier- 


Capacitor section 
from generator 


^ J Timi ng circuit 


Diode 


I Battery 


(b) Sample circuit for executing cycle A. 
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(c) Sample circuit for executing cycle B. 


ator (ref. 2) discussed in this 
report. 

The principle of operation 
of this device is best illustra- 
ted by considering a simple ver- 
sion^ such as a parallel plate 
capacitor made from a material 
whose dielectric constant varies 
with temperature. The charge on 
such a capacitor is given by the 
product of capacitance and volt- 
age. The dielectric constant can 
be decreased by either increasing 
or decreasing the temperature de- 
pending on the material employed. 
If the capacitance is decreased 
by changing the temperature of 
the dielectric while charge is 
held constant, the voltage on the 
capacitor must increase. Since 
the electric energy stored in the 
capacitor is proportional to the 
product of charge and voltage, it 
will also increase. Some of the 
thermal energy transferred in 
changing the temperature of the 
dielectric has been converted to 
electric energy. 


Figure 1. - Electrical and mechanical details of one possible design for Like any heat engine, this 

thermoelectrostatic generator. device ma y -j^g analyzed by stand- 

ard thermodynamic techniques. 

The usual expressions for internal energy and entropy must be modified to take 
into account the fact that electric and thermodynamic parameters are equally 
important in determining the performance of the generator. 


For the sake of simplicity, the two thermodynamic cycles considered are 
made up of processes in which one of the parameters of charge, electric field, 
or entropy is held constant. Cycles including isothermal processes were con- 
sidered impractical for space application because of the necessity of some sort 
of uniform rotation of the generator with respect to the radiant heat source 
(sun or reactor) in order to alternately heat and cool it. The continuously 
varying heat flux would make execution of an isothermal process difficult. 

Other possible cycles have been omitted on the grounds of lower efficiency or 
of impracticality . 


The dielectrics were selected on the basis of their potential for yielding 
the best possible performance for this generator. In order to keep the anal- 
ysis as straightforward as possible, it is desirable to find simple mathemati- 
cal approximations to the behavior of these materials. The approximations 
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used were selected on the basis of information available in the literature. 

The output parameters , thermal efficiency and energy output per unit 
volume of dielectric per cycle, were calculated as functions of maximum and 
minimum temperatures and electrical parameters. In order to arrive at the 
final parameter of interest (specific weight), it was necessary to select a 
specific design. The design selected is cylindrical (fig. l(a)) and could be 
heated by either the sun or a nuclear reactor. The heat-transfer equations 
were written for this design and solved numerically. This calculation gave the 
maximum and minimum temperatures as functions of the number of cycles executed 
per second and of the amount of dielectric material employed. The results of 
the thermodynamic and the heat-transfer calculations were combined to yield 
estimates of specific weight for this generator. 


ANALYSIS 
Basic Theory 

Energy in a variety of forms can be converted to electric energy by means 
of a variable capacitance device (all symbols are defined in appendix A). A 
parallel plate capacitor Cq is initially charged to a potential Vq so that 
the charge q on the capacitor becomes 

<1 = C-jV-l (l) 

while the electrical energy initially stored is 

w{ - | qv x (2) 


The source of the charge is then removed. The charge will remain constant 
while the capacitance is changed to a new value Cg by adding or rejecting 
heat. These conditions yield the new voltage 



and a new electrical energy 

W 2 - \ qy 2 


( 3 ) 


( 4 ) 


The net change in stored electric energy is 


W 5 


Wi 



( 5 ) 
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which represents an increase of C 2 < Cj_. Thus, some portion of the energy ex- 
pended in decreasing the capacitance is converted to electrical energy. 


For parallel plates, capacitance is given by 

( 6 ) 

which indicates that capacitance can be varied by changing either the geometry 
(for example, by forcing the plates apart and thereby doing mechanical work 
against the field) or the dielectric constant. Electrostatic generators, which 
convert mechanical energy to electric energy by varying the geometry of capaci- 
tors, have been constructed and are discussed in the literature (refs. 3 to 6) . 
If a dielectric material whose permittivity is a function of temperature is 
employed, electric energy can be produced by varying the temperature of the 
material; that is, thermal energy can be converted directly to electric energy. 
A generator operating on this principle can be analyzed by standard thermody- 
namic methods by selecting a suitable thermodynamic cycle. 



Thermodynamic Cycles 

In selecting a thermodynamic cycle, it is necessary to strive for physical 
realizability and best possible thermal efficiency. An essential requirement 
for space application is minimum specific weight; this requires the maximum 
possible energy output per unit volume of dielectric per cycle. Because of the 
necessity for minimum weight, it is desirable that the auxiliary equipment 
needed to execute the cycle be as simple as possible. Also, the output volt- 
age should be greater than the input voltage in order that a simple device, 
such as a battery, can provide the necessary excitation (i.e., initial charge 
on the capacitor) . This excitation device should not be required to deliver 
any net power over a complete cycle. 

One possible generator design is presented in figure l(a) (from ref. 2). 

It consists of a solid dielectric in the form of a thin cylindrical film with 
thin metallic electrode sections on the inner and outer surfaces. Each capaci- 
tor section is taken narrow enough that all of it may be considered to be at 
the same temperature. This device is heated by the sun and rejects heat by 
radiation; it rotates in order to carry each section through a thermal cycle. 

Execution of an isothermal process would be impractical with the rotating 
cylindrical design; consequently, the use of a Carnot cycle is not feasible. 

The remaining quantities that can be held constant during a thermodynamic 
process are charge, voltage (and therefore electric field), and entropy (assum- 
ing a lossless dielectric) . 

On the basis of the previous considerations and ease of analysis, two 
cycles have been selected for detailed analysis (see fig. 2) . The cycles are 
sketched on a graph of electric field E as a function of electric displace- 
ment D because, as will be demonstrated later in the discussion, the area 
bounded by the cycle is equal to the net electrical work output. 
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Cycle A (from ref, 2) . - An initial charge 
is placed on the capacitor at an initial field 
and an initial temperature T^ correspond- 
ing to a maximum value of the dielectric con- 
stant e. At constant charge, the temperature 
is changed to Tg so as to reduce e and, hy 
using equations (3) and (6), increase the elec- 
tric field to Eg. Next the capacitor is dis- 
charged (through the load) adiabatically (and 
isentropically for a lossless dielectric) to the 
original field E*j_. This condition results in 
a temperature T3 for which e will he slightly larger than that at Tg hut 
not as large as at T^_. Finally, when E is held constant, the temperature 
is returned to its initial value hy simultaneously adding charge and trans- 
ferring heat. 

A typical electric circuit for executing this cycle is shown in fig- 
ure l(b) . The capacitor is initially charged through the diode hy the battery 
to state 1. With both the diode and the silicon control rectifier (an elec- 
tronic switch that can he triggered on hut not off (ref. 7)) blocking current 
flow, the voltage on the capacitor is increased hy changing the temperature 
from T^ to Tg, thereby reducing e. The silicon control rectifier is then 
triggered, which allows the capacitor to discharge through the load and the 
battery (state 2 to state 3). Proper choice of load size will allow the bat- 
tery to receive exactly the energy it initially delivered. Finally, the capac- 
itor is returned to its initial state with the battery again providing the re- 
quired charge. 

Cycle B . - The first process is the same as in cycle A; that is, at con- 
stant charge the field is increased to Eg while the temperature changes to 
Tg. Electric energy is then removed (delivered to the load) adiabatically 
(and isentropically) ; this condition decreases the field to E3 and changes 
the temperature to T3. At constant charge heat is transferred, and therefore, 
the temperature is changed to T4 so that e increases, and the field drops 
to E4. Finally, the capacitor is charged adiabatically (and isentropically) 
and is returned to the initial electric field and temperature. 

A typical circuit for executing cycle B is presented in figure l(c). Dur- 
ing the first process, current is not allowed to flow from the capacitor, 
which allows voltage to be increased at constant charge. Next the bistable 
switching circuit (flip-flop) is triggered on, which allows the capacitor to 
discharge through the load and the battery; the inductor can be selected so 
that the battery receives exactly the same energy it initially delivered. When 
the field has dropped to E3, the flip-flop is triggered off, the capacitor is 
therefore isolated from the load and the battery, and the second constant 
charge process is allowed to take place. When the field has dropped to E4, 
the silicon control rectifier is triggered and the battery is allowed to charge 
the capacitor back to its initial state. 

These circuits are given as an aid to understanding the cycles and should 
not be considered the best available. Further details of electric circuitry 





(a) Cycle A. 


(b) Cycle B. 


Figure 2. - Two cycles selected for detailed 
analysis. 
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and physical construction are purposely avoided in this analysis. Another 
(less efficient) thermodynamic cycle and a variety of possible circuits can be 
found in references 8 and 9. Other possible cycles have been omitted on the 
grounds of lower efficiency or of impracticality . 


Mathematical Description of Cycles 

The processes making up the two cycles discussed can be analyzed in terms 
of the thermodynamic parameters of internal energy U and entropy S of the 
dielectric. When performing this analysis the following relations are needed: 
the first law of thermodynamics. 


dU = dQ + dW (7) 

where Q is the thermal energy; the definition of entropy, 

as = f (8) 

and an expression for dW, the electrical energy added, in terms of appropriate 
independent variables and the properties of the dielectric material. 

If the electric field E can be expressed as a univalent function of the 
electric displacement D and if the dielectric is an isotropic medium, then 
the electric energy density is given by (ref. 10) 

dW = E ID (9a) 

leading to the expression 

W = y*E(D)dD (9b) 

(all quantities are expressed in rationalized mks units) . 

When the previous expressions are used, the processes making up the two 
cycles discussed can be analyzed. In integrated form, the first lav of ther- 
modynamics becomes 


AU 


/ b y~b 

dQ + / 


E dD 


( 10 ) 


where both integrals are line integrals. 

In a constant-charge process, there will be no electrical energy trans- 
ferred; hence, the second integral in equation (10) is zero. This is true for 
a parallel plate capacitor, since the electric displacement D is numerically 
equal to the charge density. The heat transferred in this process can then be 
obtained from the end point values of U 
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dQ = U(Tb,E b ) - U(T a ,E a ) 


(Ha) 



with the condition that 


D(T b ,Eb) = D(T a ,E a ) 


(lib) 


For an isentropic process, the electrical energy transferred is given by 


f 0 E dD - UCaj^iy - U(T a ,E a ) 


(12a) 


with the condition 


- S(T a ,E a ) (12b) 

This adiabatic (and isentropic) process can be achieved in practice by 
simply charging or discharging a capacitor rapidly (assuming a lossless di- 
electric) . Such a process must result in a change in temperture for a dielec- 
tric whose permittivity varies with temperature if equation (l2b) is to be sat- 
isfied. Intuitively it would be expected that such a temperature change would 
be small. This effect has been measured (ref. ll) , for a material whose di- 
electric constant changes very rapidly with temperature, at fields up to 
0.75X10 6 volts per meter. The maximum temperature change recorded was about 
0.5° K. This small change in temperature corresponds to a very large change in 
electric field and is important in the calculation of work and efficiency; 
hence, it cannot be ignored in the analysis. 

In a constant -field process, the electrical energy transferred can be cal- 
culated directly since 


E dD * E[p(T b ,E) - D(T a ,E)] (13a) 

dQ * U(T b ,E) - U(T a ,E) - EjpC-T^E) - D(T a ,E)] (l3b) 

QThe condition demanded by the process is 

E = E a = E^ = constant 

Explicit expressions for U and S are derived in appendix B. 

Equations (ll) and (13) can be used to calculate the thermal efficiency 

heat in - heat out 

for both cycles with the usual definition, T) = heat in * This is 


(13c) 
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done explicitly in appendix C. 

When putting these processes together to form a cycle, it is necessary to 
choose a set of end point variables. Each end point is specified by a field 
and a temperature. Each process gives one equation relating two end points 
(four variables). Cycle A, which consists of three processes, is specified by 
three end points (six variables). The processes furnish three equations relat- 
ing the six variables; hence, only three variables may be chosen arbitrarily. 

By the same argument it can be shown that four variables must be chosen arbi- 
trarily to specify cycle B. 

The basic criterion used for deciding which variables to specify was ease 
of interpretation. The variables specified for cycle A were maximum tempera- 
ture, minimum temperature, and maximum electric field Eg. This choice of tem- 
peratures facilitates comparison of cycle efficiency with Carnot efficiency 
while the maximum electric field is limited by the breakdown strength of the 
particular dielectric used. For cycle B the same three variables were speci- 
fied, with the addition of the field E3 remaining in the dielectric after 
power has been delivered to the load. The ratio of this field to the maximum 
field E3/E2 gives a qualitative indication of how much energy remains in the 
capacitor after power has been delivered to the load. 


Internal Energy and Entropy 

In developing expressions for internal energy and entropy, equations ( 7 ), 
(8), and (9) are needed. All calculations that follow are subject to the as- 
sumption of a parallel plate capacitor in which edge effects are neglected. 
This assumption reduces the problem to a one -dimensional form. 

The dielectric properties of some materials, such as ferroelectrics 
(ref. 12), depend on both temperature and electric field. This relation be- 
tween dielectric constant, temperature, and field can be expressed in any one 
of the following ways : 


e = e(T,E) (l4a) 

or 

D = D[e(T),E] (14b) 

where e(T) describes the behavior of the dielectric at zero field, or 

D = e 0 E + P(E,T) (14c) 

where P is the polarization of the dielectric. For the purposes of this re- 
port, expression (lib) is used since it is a convenient form for a class of 
materials of interest. In this case T and E are taken as the independent 
variables, which gives 
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-,-p. SD de 5D 

dI) = 5FdT dT + 5E dE 


( 15 ) 


From equations (l4) and (15) it is possible to derive the following equa- 
tions for S and U (appendix B) : 


U = pcT + J' E || dE + T 


f 2* 

J IT 


3d 

Si - 


dE 


S = pc In T + 


f 2L 

J dT 


3D 


dE 


(16a) 

(16b) 


For materials such as polyethylene terepthalate, the dielectric constant 
is a function only of temperature; hence. 


and 


€ = € (T) 


D =» €e(T)E 


(17a) 


(17b) 


The functions S and U can be calculated by simply substituting these ex- 
pressions into equations (16), which yields 


r de -| 

U = pcT + [e(T) + T ^ 


q „ . e o E de 

S = pc In T + — — 


(17c) 

(I7d) 


Choice of Specific Materials 

The basic criteria for selecting specific materials for detailed study 
are the following: First, the dielectric constant must vary with temperature, 

the more rapidly the better; and second, the material must be able to store a 
large amount of electric energy per unit volume. The second criterion can be 
interpreted with the help of the expression 

w = \ e 0 £ E 2 (18) 

which is an approximation to the exact expression given in equation (9b) . This 
shows that the second criterion demands high dielectric constant and/or the 
ability to support high electric fields. 

On the basis of these criteria two materials, polyethelene terephthalate 
and barium titanate (BaTi 03 ), were selected for detailed study. The tempera- 
ture dependence of the dielectric constants of these materials is shown in 
figures 3(a) and 3(b), which were taken from references 12 to 14. (The hyster- 
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(b) Polyethylene terephthalate. 

Figure 3. - Temperature dependence of dielectric constants of two materials of 
interest for generator. 



Figure 4. - Temperature depend- 
ence of breakdown field for bar- 
ium titanate (ref. 17). Speci- 
men thickness, 0. 1 centimeter. 


esis for barium titanate shown in fig. 3(a) will be discussed later.) In addi- 
tion to the temperature variation, reference 15 shows that for barium titanate 
the dielectric constant decreases rapidly with increasing electric field. The 
maximum electric fields that can be supported are about 5x10® volts per meter 
for polyethylene terephthalate (ref. 16) and about 10^ volts per meter for 
barium titanate. The actual breakdown voltage for barium titanate is a func- 
tion of temperature (ref. 17). Figure 4 shows typical data from this refer- 
ence. The figure of 10^ volts per meter is an upper limit. 

The use of typical values of e from figure 3 (viz., 3.5 for polyethylene 
terephthalate and 10 000 for barium titanate) and the previous breakdown data 
in equation (18) indicates that both materials are capable of storing approxi- 
mately 4 joules per cubic centimeter, which is about the highest that can be 
achieved with solid dielectrics (refs. 18 and 19). When data from figure 3 are 
used in equations (5) and (6) and when both materials are assumed to give the 
same initial energy density, it can be estimated that barium titanate should 
deliver 5 to 10 times more energy per unit volume than polyethylene tereph- 
thalate for a given temperature range due to the very rapid variation of its 
dielectric constant with temperature; barium titanate would therefore lead to 
the generator with the better thermal efficiency. 


Mathematical Description of the Materials 
In order to use equations (ll), (12), (13), and (16) to analyze this gen- 
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erator, it is helpful to have a simple mathematical description of the proper- 
ties of the previous materials in the form of either equation (14a) or (14b) . 

In view of the assumptions used in the derivation of equations (l6) and (l7), 
it is necessary that the functions selected be twice differentiable; this ex- 
plicitly rules out discontinuities such as would be suggested by figure 3(a) . 
Since the purpose of this report is to find the upper limits on the capabil- 
ities of a power conversion device, as well as to show the effects of the ap- 
propriate parameters on its performance, the mathematical approximations have 
been chosen so as to produce optimistic results whenever an exact fit to the 
data is not possible. This requires that the functions selected predict a more 
rapid variation of dielectric constant with temperature and/or a greater magni- 
tude of dielectric constant than is actually shown in the data. 

Polyethylene terephthalate . - The data in figure 3(b) are very nearly 
linear over a range of temperatures; hence, the mathematical approximation 
chosen was a linear extrapolation: 

e(T) « a + (3T (19) 

where 

a = -1.6 

P = 0.0135° K" 1 

and T is expressed in °K. This provides an accurate curve fit at tempera- 
tures between 355° and 400° K. 

Barium titanate . - From the data presented in figure 3(a), it is clear 
that a simple mathematical function cannot be found that will describe barium 
titanate over the entire temperature range presented. Additional difficulties 
in finding a simple mathematical description are introduced by three properties 
of barium titanate. First, in the temperature range between 275° and 395° K, 
barium titanate exhibits a spontaneous polarization (ref. 12) . The polariza- 
tion varies with applied field E in the same manner as does magnetic induc- 
tion with magnetic intensity for soft iron (refs. 12, 13, and 20); that is, 
the variables are related by a hysteresis loop which cannot be described by a 
univalent function as required in the derivation of equations (16), (l7c), and 
(l7d) . Second, single crystals of barium titanate exhibit an abrupt change in 
dielectric constant (fig- 3(a)) at about 385° K (referred to as the "Curie 
temperature" in the literature), which corresponds to a sudden change in cry- 
stal structure (ref. 12) . The exact temperature at which this transition oc- 
curs depends on the direction of approach to the critical region (ref. 2l); 
this is indicated by arrows in figure 3(a). In addition, the Curie temperature 
is shifted upward by the application of an electric field (refs. 22 and 23). 
Finally, the properties of a single crystal of barium titanate are strongly 
anisotropic (ref. 12); the dielectric constant changes by more than an order of 
magnitude as crystal orientation is changed. 

In order to keep the calculation optimistic, an approximation based on the 
properties of a single crystal must use the properties associated with the 
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"best" orientation of the crystal; namely, the orientation that results in the 
highest dielectric constant. 

In order to avoid the spontaneous polarization and resulting hysteresis 
loops, it is desirable to limit the mathematical description to the range of 
temperatures above the Curie temperature. An analysis of a power conversion 
device using barium titanate in the region of the Curie temperature has been 
published (ref. 24) in which an attempt has been made to include the disconti- 
nuity of the dielectric constant and the dependence of the Curie temperature on 
applied field as described in reference 23. The resulting expressions, which 
involve a unit step function, are considered too unwieldy to be of use in this 
analysis. 

In polycrystalline (ceramic) form, barium titanate does not exhibit the 
discontinuity and the temperature hysteresis apparent in figure 3(a) and dis- 
cussed in reference 25. The dielectric constant is much smaller than for a 
single crystal and would result in lower performance. A theory has been pub- 
lished that predicts the behavior of such ceramics (ref. 26); however, it in- 
volves a statistical analysis that is not readily adaptable to this calcula- 
tion. 


A phenomenological theory has been proposed (ref. 27) that predicts the 
behavior of single crystals of barium titanate over the entire temperature 
range shown in figure 3(a) . This theory consists of a power series expansion 
for the thermodynamic free energy function in terms of the polarization P. 

For materials with high dielectric constant, P is approximately equal to D. 
This can be shown by examining the equation 

D » e Q E + P(E, T) - e Q e (T,E) E 

If e » e 0 , the term e Q E becomes negligible and P ~ D. In general form 
this power series is too complex to be useful in this analysis; however, it is 
shown in reference 28 that above the Curie temperature the equation can be 
simplified and rewritten in the form 


E * iyW + td3 ( 20a ) 

where 

e(T) = J _ - Tq (20b) 

D = D[e(T),E] 

if the p5 (and therefore the D 5 ) term is neglected. 

Equation 20 (a ) } which was also obtained by fitting of experimental data in 
reference 15 prior to the development of the theory of reference 28, was se- 
lected to describe barium titanate in the present report. Specific values of 
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p and T 0 (the so-called characteristic temperature) were taken from refer- 
ence 12: 


p = 170 000 
T q = 373° K 

Under the condition of low electric field, the TD 3 term of equation (20a) be- 
comes negligible and the equation can be rewritten in the form of equation 
(l7b) ; hence, equation (20b) adequately describes the dielectric in this situ- 
ation. Equation (20b) is plotted on figure 3(a) (dashed curve); it very nearly 
coincides with the data (taken at low applied field) in the region above the 
Curie temperature. For higher fields, equation (20a) must be used. Equation 
(20a) has been used in the analysis of a power conversion device (ref. 29) by 
using a linear dependence of X on temperature (as was shown to be approxi- 
mately valid in ref. 28) . The resulting expressions become unwieldy if applied 
to the thermodynamic cycles of interest in this report. According to refer- 
ence 28, T varies by about a factor of two over a temperature range of 30° K, 
a range of interest in the analysis. Published data (refs. 28 and 15) give 
typical values of T for different samples of the same material (barium tita- 
nate) ranging all the way from 5X10 9 to 5X10 1 - 1 - (rationalized mks units) . Since 
the variation in the published data is two orders of magnitude, while the tem- 
perature variation is only a factor of two, it was considered pointless to at- 
tempt to include the temperature effect. The approach taken was to consider 
T a constant and determine the effect of different (constant) values of this 
parameter . 


Some Properties of Final Equations 

The final equations describing the two cycles for both materials are ob- 
tained by substituting equations (l9) and (20) into equations (16) and (l7) , 
which are in turn used in equations (ll) to (13) . Some details of this process 
are given in appendix C. 

From the forms of some of the final equations, obtained by substituting 
the explicit expression for U and S for a given dielectric material 
(eqs. (CIO), (Cll), (C13), and (C14)), it is possible to obtain a qualitative 
idea of the performance of these cycles. An examination of these equations 
allows determination of limitations on electric field other than breakdown, of 
some relations among the endpoint variables, and of the maximum thermal effi- 
ciency for both cycles. 

Limitations on electric field . - For both cycles, the first process re- 
quires conversion of thermal energy to electrical energy at constant charge, 
which requires that e(T 2 ) < sCT-l), as indicated by equations (5) and (6) . 

This means that for polyethylene terephthalate Tg < Tpj therefore, heat must 
be rejected during the process. For barium titanate (described by eq. (20)) it 
is necessary that Tg > Tq, which requires that heat be added. Limitations on 
electric field that might be imposed by the conditions can be checked by exam- 
ining the equation 


13 


I 



r 


dQ = pc(T 2 ~ T-j_) - 


£ Eg 

o d 


JpM 

2 




< t i) + ^) | - + <§) 


'*2 


(21) 


which results from equations (ll), (lib), and (l7a) 


When equation (19) is used for polyethylene terephthalate, equation (2l) 
reduces to 


/ 


&Q = pc(Tg - T^) + 


lA [p 2 ( t i " T z) R T 1 + T 2 ) + 

2 1 (a + PT X ) 2 J 


< 0 


This sets the following upper limit on Eg : 



2pc(a, + PT-^) 2 


e o P [a^! + T 2 ) + Pp^ 



If Eg is greater than this value, the energy output will be negative; that 
is, the thermodynamic cycle will become a refrigeration cycle. Using numbers 
appropriate to polyethylene terephthalate (appendix C) gives a limit on E of 
about 3X10 9 volts per meter (T^ = 400° K, Tg = 350° K) . Since the net electri- 
cal energy generated will be zero for Eg = 0 and also for the maximum allow- 
able value of Eg (as calculated), there will be some value of Eg between 
these extremes that will yield a maximum energy output per cycle. This was 
calculated and found to be about 2x10^ volts per meter for the same T^ and 
Tg used previously. Since these values are 5 to 10 times the breakdown field 
for polyethylene terephthalate, the limitation on Eg will be the breakdown 
field rather than the field at which the energy output goes to zero. 

For barium titanate, described by equation (20b), equation (2l) reduces to 

2 

dQ = pc(T 2 - T x ) > 0 



which demands only that Tg > Tj_. There is no dependence on E (terms in- 
cluding Eg cancel out of the expression) ; hence, there is no upper limit on 
its magnitude other than breakdown. 


Relation among endpoint variables . - For 
tropic) process must result in energy output; 


both 

that 


cycles the second 


is. 



dW < 0. 


(isen- 
By argu- 
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ments similar to those used previously, it can be shown that this requirement, 
the condition on the first process, and the nature of the cycles lead to the 
following relations among the endpoint temperatures: 

Polyethlene terephthalate : 


Barium Titanate: 


Cycle A: T^_ ^ ^3 ^ ^2 

Cycle B: > T ^ > Tg 

t 4 > t 3 > t 2 


Cycle A: Tg ^ ^ ^*1 

Cycle B: Tg ^ *^1 ^ ^4 


T 2 > T 3 > t 4 


The field relations are the same for both materials: 

Cycle A: ^2 ^ = 

Cycle B: Eg > E^ > E^ 


E 2 > E, > E 4 


Maximum theoretical efficiency . - It is possible to find the maximum theo- 
retical efficiency of cycle A by using the approximate relation for barium 
titanate given by equation (20b) . Using this expression for € in equa- 
tions (ll), (12), (13), and (17) leads to the expression 


(2T n - TjT, 


T 2 - T 3 


T T, 
o 1 


n = 


T, + T_ - 2T 
13 o 


In 


T, 

c 

T, 


T 2 " T 1 


( 22 ) 


where T^ is obtained from 


e 2 _ 2pc 
2 € oP 


( T 2 - T o)' 


fa - T„\ 2 n T 3 


(23) 


1 - 


^ T 3 - 


which comes from equations (llh) , (l2d ) , and (13c) . From equation (23) it can 
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be seen that T3 = Tq gives an infinite Eg, which can be shown to give maxi- 
mum efficiency for this material. Making the substitution T 3 = Tq in equa- 
tion (22) yields 


Tn 


i T2 

In sr- 


Vax = 1 " T 2 - Tq “ Tq 


For Tg - Tq « Tq this reduces to 


1 T 2 " T 1 


j max 2 T-j 


The Carnot efficiency is given by 




T - T 
£ 1 


Carnot 


Hence , 


3 1 ^2 w 1 

T ) Carnot ^ 2 T 1 2 


‘max 


for « 1 


It is possible to find the maximum theoretical efficiency of cycle B by 
using the approximate description of barium titanate given by equation (20b) . 
The use of this expression for e in equations (17a ) y (lla) , (C14), and the 
definition of efficiency leads to the result 

i T 5 - T 4 
H *“ m rp 

2 1 

Further manipulation of equations (C14) yields the relation 


Therefore, 


(24) 


- % 


^ = ~ 


( 25 ) 


Carnot efficiency is given by the expression 


T - T 
1 2 4 


Carnot 
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since 


T„ » T . 

4 mm 

The limiting condition that must occur in order to yield maximum theoreti- 
cal efficiency for cycle B is that the net energy output approach zero. This 
equation can he written explicitly hy using equations (C2a), (C13), and (C14) 
noting that 

W = 0 = AQq^g Hh AQg ^ = AWg + AW^^j. 


This equation simplifies to the relation 

Tg - T x + T 4 - T 3 * 0 (26) 

Solution of this equation with equation (24) leads to the relations 

T 3 = T 4 (27a) 

T 2 * (27b) 


Consequently, the expression for efficiency (eq. (25)) reduces to 


il = 


T - T 
2 4 


T 2 


which is equal to the Carnot efficiency. 

The relations among the other parameters demanded by this maximum effi- 
ciency condition can be found by using equations (27 ) j, (C2a), and (C2b) with 
equations (C14) . The results of this calculation are as follows: 


% - E 4 

2gc % 4 

e oP n T 4 (T 2 - T q ) 2 " (T 4 - T 0 ) 2 

It can be seen that there is no requirement of infinite fields as was the 
case for cycle A. 

From the previous calculations it is apparent that the best possible per- 


4 


C 
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Electric field, E, V/m 


TABLE I. - TYPICAL PERFORMANCE DATA 


Figure 

Thermal 

efficiency, 

Tb 

percent 

Electric 

energy 

output, 

J/cu m 

Maximum 

electric 

field 

strength, 

^max- 7 

V/m 

Discharge 

fraction, 

e 3 /e 2 

Correction 

coefficient, 

r, 

(V)(m 5 )/C 3 

5(a) 

0.0704 

O. 0501X10“ 6 

3. 43X10 8 

— 

- 

5(b) 

.1276 

.090 

3. 43X10® 

0. 83 

- 

5(c) 

.48 

.17 

7. 22X10 6 

— 

0 

5(d) 

.78 

. 213 

7. 22X10 6 

.60 

0 




(a) Cycle A using polyethylene terephthalate. 
Maximum electric field, 3.43x10 s volts per 
meter. 


(b) Cycle B using polyethylene terephthalate. 
Maximum electric field, 3.43x10 s volts per 
meter- discharge fraction, 0.83. 



(c) Cycle A using barium titanate. Maximum (d) Cycle B using barium titanate. Maximum 
electric field, 7. 22xl0 6 volts per meter; electric field, 7. 22xl0 6 volts per meter; 

correction coefficient, 0. discharge fraction, 0.6; correction coeffi- 

cient, 0. 


Figure 5. - Typical examples of both thermodynamic cycles with both materials. 
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formance for this generator can he 



Figure 6. - Typical efficiency curves for cycle A with barium titanate. 
Maximum electric field, 7.22X10 6 volts per meter; correction coeffi- 
cient, 0. 


obtained from cycle B by using 
barium titanate as the dielectric. 


RESULTS AMD DISCUSSION 

With the preceding theory and 
approximations, a numerical anal- 
ysis was performed with an 
IBM 7094 computer. The systems of 
equations (Cll) and (C14) were 
solved by the Newton-Raphson method 
(ref. 30), and the results were 
then used in the sets of equa- 
tions (CIO) and (C13) . A thermo- 
dynamic analysis was done for 
cycles A and B by using poly- 
ethylene terephthalate and barium 
titanate. Typical numerical re- 
sults for each cycle are summarized 
in table I and figure 5 (the adia- 
batic processes are shown as dashed 
curves since their exact shape was 
not calculated) . 


The heat-transfer equations were then written for the rotating cylinder 
(fig. l(a), p. 2) in space, heated by the sun and by a nuclear reactor, and 
solved numerically (appendix D) . When the results of the thermodynamic analy- 
sis and the heat-transfer calculation were combined with the help of a suitable 
temperature scaling technique (discussed later), it was possible to determine 
the best (lowest) specific weights that could reasonably be expected for the 
dielectric material in this type of generator. 


Thermodynamic Calculations 

Cycle A . - The performance of cycle A using barium titanate (as approxi- 
mated by eq. (20a)) is shown in figures 6 to 10. The curves in figure 6 ex- 
hibit maxima which indicate that for any combination of the parameters there 
is some optimum value of T^^ - ^min "^at yields best efficiency. Figure 7 
shows the magnitude of these maxima as functions of the parameters T ^ n , 
and T, while figure 8 locates the values of T max required to yield these 
maxima. It can be seen that the best efficiencies attainable with a realistic 
value of E max (7x10® v/m) are in the neighborhood of 1.0 percent, while real- 
istic values of T m ^ n and T (393° K and 10 9 (v)(m 5 )/C 3 , respectively) reduce 
this to about 0.1 percent. The maximum theoretical efficiencies (half the 
Carnot efficiency) corresponding to these values are 1.9 and 6.4 percent, re- 
spectively. 

The energy output curves for barium titanate have the same form as the 
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Figure 7. - Maximum thermal efficiency as function of maximum electric 
field for cycle A using barium titanate. 



Figure 8. - Maximum cycle temperature required to yield maximum thermal 
efficiency for barium titanate using cycle A. 


20 




efficiency curves (fig. 6), Figures 9 and 10 give the locations and magnitudes 
of the maxima and therefore are similar to figures 7 and 8. It can he seen 
that the “best energy output available with a realistic value of (7X10 6 

V/m) is 0.45x10 s joules per cubic meter per cycle, while inclusion of realistic 
values of and T (393° K and 10 9 (V)(m a )/C 3 , respectively) reduces this 

to 0.16xl0 S joules per cubic meter per cycle. 

Cycle A was also analyzed using polyethylene terephthalate as approximated 
by equation (19). The results of this calculation can be summarized by compar- 
ing figures 5(a) and 5(c) (see p. 18) and examining the performance data 
(table I) for these cases: 

(1) Thermal efficiency, t], 0. 0704 percent; electric energy output, W, 

0.0501x10 s joules per cubic meter per cycle (fig. 5(a)) 

(2) Thermal efficiency, q, 0.48 percent; electric energy output, ¥, 

0.17x10 s joules per cubic meter per cycle (fig. 5(c)) 

As can be seen, both efficiency and energy output are considerably lower for 
polyethylene terephthalate than for barium titanate. Consequently, a more de- 
tailed discussion of this combination is not of interest. 

Circle B. - The performance of cycle B using barium titanate is presented 
in figures 11 and 12. An examination of figure (ll) shows three general fea- 
tures of interest. First, as the temperature difference T max - T m ^ n de- 
creases, the thermal efficiency increases toward Carnot efficiency. Second, 
efficiency has a parabolic dependence on the discharge fraction, the maximum 
occurring for E^/Eg - 0. Third, to a good approximation, efficiency varies 
directly as the square of the parameter Qnax* From this set of curves it can 
be seen that by using a reasonable E max (7.22X10 S V/m) , a reasonable value for 
E3/E2 (0.2), and a value of T^n (393° K) , one could expect up to 1.7-percent 
thermal efficiency (very nearly Carnot efficiency for the value of T max 
(403° K) selected). If it is assumed that x has roughly the same effect on 
this cycle as it had on cycle A (fig. 7), inclusion of a realistic value (l 0 s ) 
would reduce the expected efficiency by a factor of 5 to 0.34 percent. This 
performance is about three times better than that of cycle A. 

The energy output curves are shown in figure 12. As can be seen from 
figures 12(a) and 12(b), these curves exhibit maxima, indicating an optimum 
value for the parameter T^^, an( ^ P a ^s through zero at a value of T^ ax cor- 
responding to Carnot efficiency (fig. 11(a)). Figures 12(c) and 12(d) show 
that the net energy output varies with the parameters E^/Eg and E [nax in the 
same way as do the efficiency curves. This set of curves shows that for the 
same values of ]^ , E 3 /e 2 , and T min used in the efficiency discussion, one 

would expect an energy output as high as 0. 4x10 s joules per cubic meter per 
cycle. If it is assumed that y affects this cycle as it did cycle A (fig. 9) , 
inclusion of a realistic value (10 s ) would reduce the expected energy output by 
a factor of 2 to 0. 2x10 s joules per cubic meter per cycle. 

Cycle B was also analyzed using polyethylene terephthalate as approximated 
by equation (19). The results of this calculation can be summarized by com- 
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Thermal efficiency, tj, percent 



380 400 420 440 460 480 500 0 . 2 . 4 . 6 . 8 1.0 2 4 6 8 10 20xl0 6 

Maximum cycle temperature, T max , °K Discharge fraction, E 3 /E 2 Maximum electric field, E max , V/m 


(a) Maximum electric field, 7. 22xl0 6 volts per meter; correction (b) Minimum cycle temperature, 393° K ; maximum cycle (c) Minimum cycle temperature, 393° K. 

coefficient, 0. temperature, 408° K; correction coefficient, 0. 

Figure 11. - Efficiency curves for cycle B with barium titanate. 





Energy output per cycle, W, J/cu m 



(a) Minimum cycle temperature, 393° K; maximum electric field, 
7.22xl0 6 volts per meter; correction coefficient, 0. 



390 410 430 450 470 490 510 

Maximum cycle temperature, T max , °K 


(b) Maximum electric field, 7. 22xl0 6 volts per meter; discharge 
fraction, 0 ; correction coefficient, 0 . 



0 .2 .4 .6 .8 1.0 


Discharge fraction, E 3 /E 2 

(c) Minimum cycle temperature, 393° K; maximum cycle tempera- 
ture, 418° K; correction coefficient, 0. 



(d) Minimum cycle temperature, 393° K; correction 
coefficient, 0 . 


Figure 12. - Energy output, curves for cycle B with barium titanate. 
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paring figures 5(*b) and 5(d) and examining the performance data (table I) for 
these cases: 

(1) Thermal efficiency, r\, 0.1276 percent; electric energy output, ¥, 

0. 090X10 6 joules per cubic meter per cycle (fig. 5(h)) 

(2) Thermal efficiency, q, 0. 78 percent; electric energy output, W, 

0. 213x1 0^ joules per cubic meter per cycle (fig. 5(d)) 

As can be seen, both efficiency and energy output are considerably lover for 
polyethylene terephthalate than for barium titanate. Consequently, a more de- 
tailed discussion of this combination is not of interest. 

Best performance . - From the preceding discussion and figures it is clear 
that barium titanate is the more promising of the tvo dielectrics considered. 
When this material is used, cycle B is tvo to three times as efficient as 
cycle A and yields roughly the same energy output per cycle. A comparison of 
figures 5(c) and 5(d) and the attached performance data provides a numerical 
example of this trend. Since the goal of this study is to achieve maximum 
pover output vith minimum veight, the combination of barium titanate and 
cycle B is the only one considered in subsequent specific veight calculations. 


Heat-Transfer Calculations 


In order to find the specific veight of this system, it is necessary to 
find the number of cycles executed per second and the corresponding maximum and 
minimum temperatures, vhich in turn give the energy output per unit volume per 
cycle from the thermodynamic analysis. Specific veight is then given explic- 
itly by the expression p/Wf. This is accomplished by solving the heat- 
transfer equations (appendix D) for the design of figure l(a) and using the 
results (maximum and minimum temperatures) vith the thermodynamic data (viz., 
the energy output per unit volume curves) . 


Typical temperature distributions on a cylinder rotating in space close to 
the Earth and heated by the sun vere obtained from equation (Dl) and are pre- 
sented in figure 13. The "rotation parameter" A is given by 


A = 


ab r 


2jcfLpc 


(29) 


as in reference 2. Taking the special case of surface emissivity (Bg = 1) vith 
p and c taken at values appropriate to barium titanate (appendix C) and ex- 
pressing L in meters lead to 


A = 3. 95xlO -11 (fTi) -1 


The effect on the temperature distribution of executing cycle A vith a 
hypothetical material having dielectric properties similar to those of barium 
titanate (scaled to fit the temperature range encountered) vas calculated for 
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Figure 13. - Temperature distribution for thin -walled hollow cylinder in vicinity of Earth. Incident 
power, 1390 watts per square meter. 


a range of the parameters A and ^min “ T 0 (appendix D) . The solution shown 
in figure 14 represents a case in which the effect on the temperature distri- 
bution is a maximum. Since all thermodynamic calculations were done with 
T max and T m j_ n as parameters, figure 14 shows that the effects of the thermo- 
dynamic cycle on the temperature distribution may be neglected and the results 
for the "inert cylinder” used. This conclusion is also valid for cycle B, be- 
cause the only effect of cycle B would be to introduce a second discontinuity 
into the temperature distribution (from 1^ n to T min + 2° or 3° K) . 

Since lighter weight can result from higher operating temperatures, the 
possibility of heating the generator with a nuclear reactor was considered. 

When figure 18 and the equations in appendix D were used, temperature distri- 
butions were found for the case of reactor heating. 

The maximum and minimum temperatures occurring on the cylinder are shown 
as functions of the rotation parameter A for both reactor heating and solar 
heating of the "inert” (dW - 0) cylinder in figure 15. These temperatures are 
not, in general, compatible with the two materials used in the thermodynamic 
analysis (fig. 3, p. 10) . 


Specific Weight 

Tempe rature s caling . - There are several materials known which behave like 
barium titanate but in a different temperature range (refs. 12 and 31 to 34); 
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Figure 14. - Temperature distribution for thin-walled hollow cylinder in vicinity of Earth with effect of 
thermodynamic cycle included. Rotational parameter, 1. OxlO' 9 radian -1 °K - 3 ; cycle A; material, 
ferroelectric with characteristic temperature of 273. 25° K; correction coefficient, 0; maximum 
electric field strength, 7.22xl0 6 volts per meter; incident radiant energy from sun, 1390 watts per 
square meter. 


400 



100 1 — 



Figure 15. - Maximum and minimum temperatures on thin-walled hollow cylinder heated by reactor only and by 
both reactor and sun as functions of rotation parameter. Incident power from reactor, 300 000 watts per 
square meter; incident power from sun, 1390 watts per square meter. 
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these materials represent a range in T Q (as in eq. 20(h)) from 100° to about 
750° K. In order to use the data of figure 15 to obtain estimates of the spe- 
cific weight of a generator, it is necessary to construct mathematical models 
for materials that exhibit dielectric properties exactly like those of barium 
titanate in the temperature ranges calculated. Since T^q^. and T^j^ have 
been used as parameters in the thermodynamic analysis, it is necessary to pre- 
serve the difference T max - T m ^ n and the corresponding change in the dielec- 
tric constant. This can be done by simply "sliding” the temperature scale of 
figure 3 (p. 10) up or down to fit the data of figure 15. 

For barium titanate, described by equation (20), the tern 

e (T) = T Z t o 

determines the temperature dependence of the dielectric constant; the charac- 
teristic temperature T 0 is a convenient reference. The required scaling can 
be accomplished by simply shifting T 0 up or down to fit the desired tempera- 
ture range and noting that the value of T 0 = 373° K was used for figures 6 
to 12. Values of T - T Q less than 10° yield values of e that are beyond 
the maximum values observed experimentally. Consequently, values of T m j_ n - T 0 
less the 10° are not considered in subsequent calculations. In order to eval- 
uate the effect of this temperature scaling, cycle A was analyzed for a range 
of values of T Q using equations (C2a) , (CIO), (Cll), and (20b). The results 
for the material similar to barium titanate are presented in figure 16. From 
figure 6 (p. 19) it can be seen that ^max ~ Tftiin ~ 20° K yields maximum ther- 
mal efficiency for most cases; figure 16 shows that for this value of 
T max - ^mln ^•* le magnitude of the thermal efficiency does not change by more 
than ±25 percent over the range of temperatures of interest. 



Figure 16. - Effect of absolute temperature on thermal efficiency of cycle A 
using ferroelectrics similar to barium titanate. Correction coefficient, 0. 
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The effect of scaling the temperatures, as shown 
in figure 16, is small enough that it can he ignored 
in estimating specific weights. This approximation 
is of the same order of magnitude as other approxima- 
tions made thus far in setting up the analysis. 
Therefore, figures 12(a) and 12(h) can he used di- 
rectly in a specific weight calculation hy simply 
scaling the endpoint temperatures, as discussed pre- 
viously, to fit the data of figure 15. 


Values of specific weight . - When the tempera- 
ture scaling technique discussed previously, the data 
of figures 12(a) , 12(h), and (15) , and the definition 
of the rotation parameter A (eq. (29)) are used, spe- 
cific weights may he calculated for various combina- 
tions of rotational frequency and thickness of the 
dielectric. Figure 17 shows typical specific weight 
curves for the material similar to barium titanate 
using cycle B with solar heating. The dielectric 
thickness chosen (l. 27x10"^ m) is about the lower 
limit that can he achieved with barium titanate with- 
out a large decrease in dielectric constant (refs. 35 
and 36). The equations for A and for specific 
weight can he combined to yield specific weight = 2rrcp 2 A L/crSgW. From this form 
it can he seen that specific weight is directly proportional to dielectric 
thickness. The lowest specific weights attainable with solar heating are in 
the neighborhood of 50 to 150 pounds per kilowatt, a value comparable to that 
of solar cells (ref. 37). 


Figure 17. - Typical specific weight 
curves for cycle B with ferro- 
electric. Maximum electric field, 
7. 22xl0 6 volts per meter; correc- 
tion coefficient, 0. 


Figure 15 indicates that reactor heating would decrease specific weight of 
the dielectric material hy roughly two orders of magnitude for a given thick- 
ness, since specific weight is directly proportional to the A required to 
yield a given temperature range and hence specify a value of W. It is nec- 
essary, however, to consider the weight of the radiator required to transfer 
heat from the reactor to the generator. 


An estimate of radiator weight for the heat flux used in figure 15 was 
obtained hy using a radiator design computer program (unpublished data obtained 
from Arthur V. Saule, NASA, Lewis Research Center). This gave a lower limit of 
roughly 0. 25 pound per radiated kilowatt. Assuming a thermal efficiency of 
0. 75 percent for the generator (optimistic) leads to a minimum radiator weight 
of 33 pounds per kilowatt output. In this estimate, reactor weight and struc- 
tural materials needed for the generator have been ignored. 


CONCLUDING REMARKS 

The performance of a thermoelectrostatic generator has been analyzed in 
detail. Two thermodynamic cycles have been selected, and mathematical expres- 
sions for internal energy and entropy derived. Two materials were selected 
and their properties approximated by simple mathematical functions. The equa- 
tions resulting from this mathematical model were solved numerically. 
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The heat-transfer equations -were written for a typical design, heated *by 
the sun and hy a nuclear reactor, and solved numerically. The results of this 
calculation were combined with the results of the thermodynamic analysis to 
yield the specific weight of the system. 

From typical illustrative calculations the following features of generator 
performance were evident: 

1. Barium titanate and materials similar to it are capable of producing 
three to four times the energy output per unit volume at three to four times 
the thermal efficiency possible with polyethylene terephthalate. 

2. Thermodynamic cycle B is capable of two to three times the energy out- 
put at two to three times the thermal efficiency of cycle A. 

3. The best possible thermodynamic performance for this generator is a 
thermal efficiency of 1 to 4 percent and an energy output of 0.5 to 1.5 joules 
per cubic centimeter of dielectric per cycle, while realistic values are prob- 
ably lower by at least a factor of five. 

4. With solar heating, specific weight of the generator is limited by 
minimum allowable dielectric thickness and thermodynamic performance of the 
dielectric to about 50 to 150 pounds per kilowatt. 

5. With reactor heating, specific weight is limited by the weight of the 
radiator needed to transfer heat from the reactor to the generator and by the 
thermodynamic performance of the dielectric to about 30 to 60 pounds per kilo- 
watt. 


6. With currently known materials, this generator is too heavy to be of 
interest for electric propulsion. As a source of auxiliary power it is com- 
parable, at best, to solar cells in weight. If materials with substantially 
better characteristics than those considered are found, competitive performance 
as a power source might be possible. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, January 14, 1965. 
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APPENDIX A 


SYMBOLS 

area, sq m 

generalized endpoints of thermodynamic process 

constant defined on p. 36 

capacitance, F 

specific heat, j/(kg)(°K) 

electric displacement, c/sq m 

electric field strength, v/m 

discharge fraction 

rotational frequency of cylindrical generator, cycles/sec 

incident radiant energy from sun, 1390 W/sq m 

incident radiant energy from reactor, 300 000 W/sq m 

thickness of dielectric film, m 

polarization of dielectric, c/sq m 

"Curie constant” used in eq. (20b), °K 

thermal energy (heat) per unit volume, j/cu m 

electric charge, C 

entropy per unit volume, j/(°K)(cu m) 

entropy per unit volume at zero electric field, j/(°K)(cu m) 
absolute temperature, °K 

characteristic temperature for barium titanate, °K 
internal energy per unit volume, j/cu m 

internal energy per unit volume at zero electric field, j/cu 
voltage , V 



¥ 


¥’ 


x,y,z 


a 

P 

r 

& 2 

e 


n 

e 


a 

A 


P 

Subscripts: 


electric energy output per unit volume per cycle, j/cu m 

electric energy (total) , J 

arbitrary variables 

constant used in eq. (19) 

constant used in eq. (19), l/°K 

"correction coefficient" used in eq. (20a), (v) (m 5 )/c 3 
surface absorptivity 
surface emissivity 
dielectric constant 

permittivity of free space, 8, 85xl0“ 12 F/m 

thermal efficiency, percent 

angle 

Stefan- Boltzmann radiation constant, 5. 7x1 CT 8 w/(°K 4: ) (m 2 ) 
"rotation parameter” in heat-transfer equations, l/(°K 3 ) (rad^) 
density, kg/cu m 


a,h generalized endpoints of thermodynamic process 

max maximum value of variable occurring in course of thermodynamic 

cycle 

min minimum value of variable occurring in course of thermodynamic 

cycle 

1,2,3, etc. value of variable at endpoints of thermodynamic cycle} numbers 

refer to order in which endpoints occur in test discussion 
(i.e., both cycles begin and end at point 1) 

1-2, 2-3, etc. refer to thermodynamic processes from one endpoint of cycle to 
next 
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APPENDIX B 


DERIVATION OF THERMODYNAMIC FUNCTIONS 

In the following discussion mathematical details of the derivation of the 
expressions (l6a) and (l6b) for U and S in terms of the independent vari- 
ables T and E are presented. The basic equation needed is the first law of 
thermodynamics : 


dU = dQ + dW 


Making the substitution 


dW =* E dD 


gives 


dU = dQ + E dD 

For the particular form of D in equation (l4b). 


3D d€ 


3D 


“ " 3F df dT + st 


(Bl) 


(B2) 


When internal energy is expressed in the form U = U[e(T),E], equation (Bl) can 
be written 

dU - dQ + u(j| § dT + || as) - g § dT + | dE (B3) 

Similarly, taking S = S[e(T),E] leads to 


dQ _ .o _ 3S de 
~¥ ~ dS " 5e dT 


dT + 



(B4) 


Since dS is an exact differential, it is necessary that 


3 /3s de\ 3 /3s \ 

3e dT ) 3t ^3e J 

Solving (B3) for dQ, substituting into (B4), equating coefficients, and 
applying (B5) yield 


3 ] 

/3u de 

v d£ \ 

3 [*] 

/3u 

v Sd\ 

3e H 

p \3T dT 

E 3? dT/ 

= 3t g 

P \oE ~ 

E 3eJ 


(B5) 


(B6) 


Carrying out this operation and noting that 

3 2 u _ 3 2 u 
3e 3e = 3e 3E 
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5 2 d _ £> 2 d 

c^E 5e be 


yield 


Integration with respect to 

U = U 0 (T) + f e|» + I / i sf « 
which is equation (l6a) . 

To calculate S, it is noted that (B5) and (B6) yield 

dS _ 1 /du 5 d\ 

51 " T y dE " E 51/ 

and 

bS _ 1 /bU c>D de\ 

5f " T dT " h 5e dTj 

Substituting (B7) into (B9) gives 


c>U _ ciD de , _ bV 

“ T 5FdT + E 5 e 


OU 

51 * 


E gives 


5S _ SD de_ 
51 “ 5e dT 

so that integration with respect to E gives 


S = S 0 (T) + y g d£ 

To find 80(1), it is noted that the condition E = 0 reduces (BIO) to 

asp _ 1 d.u 0 

dT ~ T dT 


r /1 du o\ 

S o( T ) = J dT - 

For a solid material, if it is assumed that p and c are independent 
temperature, Uq(T) » pcT, and hence, 

S q (T) = pc In T 

Combining equations (Bll) and (B13) yields equation (16b) . 


(B7) 


(B8) 


(B9) 

(BIO) 


(Bll) 


(B12) 

of 

(B13) 
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appendix: c 


DEVELOPMENT OF FINAL EQUATIONS 

In this appendix the explicit equations used in the thermodynamic analysis 
are derived. The exact equations for U and S are obtained for the two di- 
electrics considered. The use of these expressions in analyzing the two ther- 
modynamic cycles is then discussed in detail. 


Internal Energy and Entropy 

For polyethylene terephthalate the expressions for internal energy and 
entropy are obtained by substituting equation (19) into equations (17c) and 
(17 d) giving 


€ E 2 


U(T,E) « pcT + - 2 — ( a + 2pT) 

(Cla) 

e E 2 


S(T,E) = pc In T + p 

(Clb) 


For barium titanate, the special case of T = 0 is considered first. 
Equation (20a) may then be rewritten 

D = e o e(0?)E 

which is the same form as equation (l7b) . Consequently, the expressions for U 
and S can be derived by substituting equation (20b) into equations (17c) and 
(I7d): 


e/Er pT 

U(T,E) = pcT - 2 __ (C2a) 

2 ( T _ Tq )2 


€ E 2 

S(T,E) B pc In T - - 2 — E — (C2b) 

2 (t - T 0 ) 2 


In the general case, T / 0 and S and U are found by substituting 
equation (20a) into equations (l6a) and (l6b) . To do this, it is convenient to 
solve equation (20a) for D (ref. 38): 
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where 



To complete the calculation of S and U from equations (l6a) and (l6b) for 
X 0, it is necessary to compute the integrals 

and 

f E 5 §“ 


Equation (C3) yields 



in equation (C4), the integral reduces to the form 



(C5) 


When the substitutions X = BY in the first integral and X = -BZ in the 
second are made,, equation (C5) reduces to a standard form in reference 39. 
This leads to the result 
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To find J* E ^ dE, it is noted that an integration by parts gives 

J E || dE = ED - y D dE 


(C6) 

(C7) 


If equation (C3) and the same substitutions as before are used, the remaining 
integral takes the form 


/x^ + B 2 + x) dX 



Evaluating this equation (ref. 39) and using equations (C3) and (C7) lead to 



tions (l6a) and (l6b) yield 
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S(T,E) = pc In T - 



[VSHSJ ♦ kf 



(C9b) 


Defining Equations for Cycle A 

In the following discussion, the endpoints of the cycle are labelled 


and 3-1 (constant charge, adiabatic, constant field), respectively, 
tions (ll) to (l3) are used, cycle A may be analyzed as follows: 

When equa- 


AQq_ 2 = U(T 2 ,E 2 ) - UCT-^Rl) 

(ClOa) 


AW 2 _ 3 = U(T 3 ,E 3 ) - U(T 2 ,E2) 

(ClOb) 


aq 3 _i « u^Ep - u(t 3 ,e 3 ) - e^dCTj^) - d(t 3 ,e 3 )] 

(ClOc) 


AW 3-1 - E i[ D (Ti,Sl) - D(T 3 ,E 3 )] 

(ClOd) 

where 

d(t 2 ,e 2 ) * dCt-l,^) 

(Clla) 


s(t 3 ,e 3 ) = s(t 2 ,e 2 ) 

(Cllb) 


E 1 = E 3 

(Cllc) 

Equations 

(Oil) yield E^, E 3 , and T 3 since T^, T 2 , and Eg are 

specified. 


In order to calculate the thermal efficiency, the usual definition 

n = k ea/ k ^- n - Q^t use( ^ For a ma terial whose dielectric constant de- 

1 heat m 

creases with increasing temperature, heat must be added during the first (con- 
stant charge) process in order to have Eg > Ej_ as required. Therefore, 

AQi _2 = heat in and AQj_j_ = heat rejected, which leads to the expression 


u(t 2 ,e 2 ) - u(t 3 ,e 3 ) - e 1 [d(t 1 ,e 1 ) - d(t 3 ,e 3 )] 

U(T 2 ,E 2 j - U(T 1 ,E 1 ) — 


(C12a) 


For a material whose dielectric constant increases with increasing tem- 
perature, heat must be rejected during the constant-charge process; hence. 
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AQi _2 = heat rejected and AQ 3-1 * heat added. This leads to the expression 


T) 


u(t 2 ,e 2 ) - u(t 3 ,e 3 ) - e ^(t-^e-l) - d(t 3 ,e 3 )] 

U^E^ - U(T 3 ,E 3 ) - E 1 [D(T 1 ,E 1 } - D(T 3 ,E 3 J] 


(C12b) 


Cycle B 

In the following discussion, the endpoints of the cycle will he labelled 
1, 2, 3, and 4, respectively, and the processes joining them will be labelled 
1-2, 2-3, etc. (constant charge, adiabatic, constant charge, adiabatic, respec- 
tively). Equations (ll) to (l3) lead to the following equations: 


a %_2 = UCT^Eg) - U(T 1 ,E 1 ) (C13a) 
AW 2 _ 3 = U(T 3 ,E 3 ) - uCTg^Eg) (C13b) 
AQ 3 _ 4 = U(T 4 ,E 4 ) - U(T 3 ,E 3 ) (C13c) 
AW^ = U(T 1 ,E 1 ) - U(T 4 ,E 4 ) (C13d) 

where 

D(T 1 ,E 1 ] « D(T 2 ,Eg) (C14a) 
S(T 3 ,E 3 ) = S(T 2 ,E 2 ) (C14b) 
D(T 4 ,E 4 ) = D(T 3 ,E 3 ) (C14c) 
S(T 1 ,E 1 ) =* S(T 4 ,E 4 ) (C14d) 


If equation (l 6 a) and the equation for D appropriate to the dielectric being 
considered are used, equations (C14) can be solved for Tj_, T^, Ej, and E 4 
since Tg, T^, Eg, and E^ are specified. 

Thermal efficiency can be calculated from equations (C13) . The arguments 
used in the efficiency calculation for cycle A apply equally well here. There- 
fore, the use of a material whose dielectric constant decreases with increasing 
temperature requires that AQ-^g =* heat added and a heat rejected; 

hence, 


u(t 2 ,e 2 ) - uCt-^eP + u(t 4 ,e 4 ) - u(t 3 ,e 3 ) 
u(t 2 ,E2) - UCT-^E-l) 


(C15a) 


If a material whose dielectric constant increases with temperature is 
used, ** heat rejected and AQ^^ - heat added; hence. 
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U(T 2 ,E2) - Ud-L,^) + U(T 4 ,E 4 ) - U(T 3 ,E 3 ) 

u(t 4 ,e 4 ) - u(t 3 ,e 3 )' 


Table II gives the densities and specific heats used for the analysis. 


TABLE II. - DENSITIES AHD 
SPECIFIC HEATS OF TEST 


MATERIALS 


Density,, 

P 

Specific heat, 
c 

kg/m 3 

g/cm 3 

j/(kg) (°K) 

cal/g 


Barium titanate 


a 5500 

a 5.5 

a 418.6 

a 0.1 

Polyethylene terephthalate 

'bisso 

^l. 38 

c 1300 

c 0.31 


^Reference 24. 

reference 14. 
C Reference 40. 


(C15b) 
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APPENDIX D 


HEAT-TRANSFER CALCULATIONS 

In this section, the equations used to describe the energy balance for the 
cylindrical generator design are presented. Two heat sources, the sun and a 
nuclear reactor, are discussed. Finally, the effect of the transfer of elec- 
trical energy is included for the special case of cycle A with barium titanate 
as the dielectric. 


Cylinder - Solar Heating 

In setting up the heat-transfer equation for a hollow cylinder heated by 
solar radiation (fig. l(a), p. 2), the following approximations are used: 

(1) Heat conduction along the dielectric film is ignored; radiation is the 
only heat-transfer mechanism considered. 

(2) No heat is transferred to the interior of the cylinder by radiation; 
the inner electrode is assumed a perfect reflector. 

( 3 ) Only thermal properties of the dielectric are considered. 

(4) Uniform rotational frequency is assumed. 

The first two of these assumptions are discussed in detail in reference 2. 

These four approximations lead to the following description of the heat- 
transfer process (ref. 2) : 

3 K 

H = A sin 0 " ttr 4 0 < 0 Jt (Dla) 

2 

H=-AT 4 rt < 0 < (Dlb) 


where 

7\ = abg/Sj tfLpC, Kj_ * 1390 watts per square meter (solar constant at Earth*s 
orbit), and 0 < 0 < jt corresponds to the sunlit side. In order to solve 
these equations, it is necessary to choose an initial temperature at an initial 
angle. The solution to one equation provides boundary conditions for the next 
as the "critical" angles 0=0 and 0 = it are encountered. 


Cylinder - Reactor Heating 

Figure 18 shows the model used to determine the effect of heating the 
cylinder with a reactor. Heat is transferred from the reactor to the cylinder 
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Radiator by the radiator. In the interest of mini- 

mum weight, the radiator was assumed to 
yield a heat flux density of 300 000 watts 
per square meter, which corresponds to a 
blackbody temperature of 1520° K. This 
temperature is a practical upper limit for 
materials presently under consideration for 
use in a radiator (ref. 4l) . The maximum 
temperatures of the dielectric compatible 
with known ferroelectric materials is in 
the range of 700° to 800° K (ref. 34). The desired heat flux density and the 
desired temperature range of the dielectric dictate that the radiator subtend 
an arc of approximately 20°, as shown in figure 18. 


Solar 

radiation 


-Generator (dielectric film) 

Figure 18. - Model used to determine effect of heating a 
cylinder with a reactor. 


When subject to the same assumptions used for solar testing, this system 
can be described as follows: 


dT ^1 ^1 

ae - * - sin 0 " 


0 < e < 80°; 100° < e < 180 c 


(D2a) 


dT . & 1 K 2 . fl 
T7T = A r- — sm e 


<3.0 


5 2 o 


80° < 0 100° (D2t>) 

(for T « 1520° K = T radiator) 


— = - At 4 

ae A 


180 < 0 < 36 0 L 


(D2c) 


Effect of Thermodynamic Cycle 

In order to determine the change in temperature distribution on the cylin- 
der (given by eqs. (Dla) and (Dlb)) caused by the execution of a thermodynamic 
cycle, it is necessary to modify the mathematical description of the dielectric 
since the temperature range on the cylinder will not in general match that in 
figure 3 (p. 10) (viz, 350° to 400° K for polyethylene terephthalate and 385° 
to 450° K, above Curie temperature, for barium titanate) . This was done by 
considering a hypothetical material with the properties of barium titanate in a 
different temperature range. The temperature scaling technique described in 
the text was used. 

In the interest of simplicity, cycle A was used in this analysis with a 
dielectric having properties similar to those of barium titanate (temperatures 
scaled). Equations (C2a), (C2b), (CIO), and (Cll) describe the thermodynamic 
cycle. 

Heat is transferred by radiation; hence, an equation of the form 

^^thermodynamic 5 * ^-^radiated mus ' t be satisfied. Because of this additional con- 
straint on the temperature, only two variables can be specified arbitrarily. 
These were E max and A. An additional stipulation was that T^ ^ - T Q = 10° K 
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so as to guarantee a maximum variation of dielectric constant and hence a maxi- 
mum effect on the temperature distribution of the cylinder. Writing equa- 
tion (ClOa) in differential form leads to the equation 


dT _ , ^1 

ae 6 2 a 


sin 0 - AT 4 


@1 9 < 


(D3a) 


(which describes process 1 - 2 ) where 0 -j_ is the angle at which T m j_ n occurs 
and 0 g is the angle at which T max occurs. 

If an adiabatic (and isentropic) process (process 2-3) is to be achieved 
in a practical system, it must be done quickly, that is, while the generator 
rotates through a small angle. Therefore, it was arbitrarily assumed that the 
generator rotates from 63 to ^2 + 2° during this process. (The calculation 
was done for several choices of this angle; results were very nearly identical 
for angles ranging from 0° to 5°.) The temperature T 3 , occurring at 63 , is 
calculated from equation (23) : 



/ 1 \ 2 


1 

O 

DO | 

2pc 

V t 2 - V 

1 - 

- V J 


+ In jp „ 0 

l 2 


(D3b) 


Process 3-1 (constant field) is described by putting equation (ClOc) in differ- 
ential form, which gives 


1 + 


o 2 
pc 


T 


T 


IT. 


2" T o / (T-Tj- 


~ = A % ^ sin 0 - AT 4 

d0 &2 0 


= - AT 4 


jt < 0 < 2 tc 



(D3c) 


Equations (D 3 ) were solved numerically. It is necessary to select spe- 
cific values for 0j_ and 0 2 , at which angles Tj_ = T^-j n and T 2 * T max must 
occur, in order that the equations can be solved. These arbitrary initial 
choices will not, in general, be correct (i.e., T mir) and T max will occur at 
angles other than those selected) ; hence, an iteration process is needed to 
find the correct values of 0j_ and 0 2 , at which T^n and T max occur, and 
to readjust the value of T q to guarantee that 


T . 
min 


T ■ 10° K 
o 


as originally specified. Equations (CIO) and (C 12 a) then yield thermal effi- 
ciency, net energy output per cycle, etc. 
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